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CONVERGENCE OF REGULARIZED NONLOCAL INTERACTION 

ENERGIES 

KATY CRAIG AND IHSAN TOPALOGLU 


Abstract. Inspired by numerical studies of the aggregation equation, we study the effect of 
regularization on nonlocal interaction energies. We consider energies defined via a repulsive- 
attractive interaction kernel, regularized by convolution with a mollifier. We prove that, 
with respect to the 2-Wasserstein metric, the regularized energies P-converge to the unreg¬ 
ularized energy and minimizers converge to minimizers. We then apply our results to prove 
T-convergence of the gradient flows, when restricted to the space of measures with bounded 
density. 


1. Introduction 

We consider the nonlocal interaction energy 
(1.1) E(n) = f f K{x-y)dn(x)dii{y), 

JR d jR d 

over the space T^O^) of probability measures with finite second moment, where the pairwise 
interaction kernel K : -a- Muj+oo} is an even, locally integrable, and lower semicontinuous 

function. Depending on the choice of K, the asymptotic states of many physical and biological 
systems can be characterized as minimizers of this energy. Of particular interest are interaction 
kernels which are repulsive at short distances and attractive at long distances. Important 
examples of such repulsive-attractive kernels are Morse-type potentials 

K(x) = C r e~\ xUlr - C a c~W /la 

(L2) _ . 

with C r /C a < ( l r /l a ) , 0 < l r < l a , and 0 <C a < C r , 


and potentials in the power-law form 

(1.3) K(x) = \x\ q /q — \x\ p /p with — d < p < q, 



Due to the competition between the repulsive and attractive terms, the minimization of 


these energies leads to complex equilibrium configurations (4||ll fl4l| 16,[32],|34^36|39,[51|, f57[[58] . 
In the case of repulsive-attractive power law kernels B local minimizers exhibit a variety 
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of qualitatively different patterns, from solid rings to broken, rounded triangles 40 , and the 


dimension of the minimizers’ support can be characterized in terms of the strength of the 
repulsive forces [6]. 

In addition to qualitative properties of local minimizers, there has also been significant inter¬ 


est in the existence and uniqueness of global minimizers. (Since (1.1) is translation invariant, 
minimizers are only unique up to translation.) Existence was recently established for a range 
of kernels, including (1.3), by Simione, Slepcev, and the second author J54j and by Cahizo, 

On the other hand, characterizing which kernels have a unique 


Carrillo, and Patacchini [18] 
global minimizer is essentially open, due to lack of convexity. In the particular case of Coulomb 
repulsion and quadratic attraction (when the energy is convex with respect to H -1 ) the unique 
global minimizer is the indicator function on a ball (c.f. for example [21,28]). 

Alongside the static problem of characterizing the minimizers of the interaction energy, 
there has also been significant interest in the dynamic problem of understanding the behavior 
of systems as they evolve toward a local minimizer. Such systems arise in biological swarming 
robotic swarming [29,49], and assembly of viral capsid proteins 38 , and they may 


47 18 58 


be modeled as gradient flows of the energy with respect to the 2- Wasserstein metric. Formally, 
this corresponds to the nonlinear, nonlocal partial differential equation 


(1.4) 


pt + V • (vp) = 0 
P(x, 0) = po(x), 


with v = —2VAT * p, 


known as the aggregation equation. (We choose to put a factor of two in the velocity field 
instead of putting a factor of one half on the energy (1.1).) For semi-convex interaction kernels 


K, with up to a Lipschitz singularity, weak measures solutions to (1.4) exist for all time and are 
unique 0H) . If the kernels are sufficiently convex, solutions converge exponentially to a unique 
global minimizer 24 . However, for nonconvex kernels with merely integrable singularities, 


much less is known about the evolution of measure solutions and their asymptotic behavior. 
In the particular case of Coulomb repulsion and quadratic attraction, solutions with bounded, 
continuous initial data converge to the unique global minimizer algebraically in time |14| . 

Due to the analytical difficulties repulsive-attractive kernels present, both theoretical and 
applied work is often complemented by numerical simulations. The most common method for 


simulating (1.4) is a particle method, in which one approximates the initial data by a sum 
of Dirac masses, po ~ where Xi is the location of the Dirac mass and ra,; ^ 0 is 


its weight. The corresponding solutions of (1.4) are formally of the form p ~ 

where the trajectories of the Dirac masses X;(tj satisfy the following finite system of ODEs: 


(1.5) 


j t Xi{t) = -2 ]T VK(W(f) - Xfltflmu 

3 = 1 


N 


Xi{ 0) = Xi. 


For a range of interaction kernels, including repulsive-attractive power-laws (1.3) in the pa¬ 
rameter regime 2 — d < p ^ 2, q > 0, Carrillo, Choi, and Hauray prove the convergence 


of the particle method solution ^ Xi{t) m i to weak measure solutions of the aggregation 


equation 20 


In recent work, Bertozzi and the first author consider a modification of these types of particle 
methods [3l], analogous to classical vortex blob methods from fluid dynamics (c.f. |[3][7j|8] and 
references therein). Specifically, since the gradient of the kernel VAT may be singular at the 
origin, they regularize the interaction kernel by convolution with a modifier, leading to an 
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approximate solution of the form YliLi where Xf(t) satisfies (1.5) with VK replaced 

by X7(K *ip e ). By regularizing the interaction kernel K, the authors extend particle methods to 
a wide range of interaction kernels, including the power-law kernels (1.3) for 2 — d ^ p q, and 
obtain quantitative rates of convergence to classical solutions. While this convergence result 
is limited to bounded time intervals, numerical results indicate that regularizing the energy in 
this way may also be useful in studying asymptotic behavior (see Figure [T]). 



initial data 




Figure 1. Numerical solutions of (1.4) at t = 100 for Coulomb repulsion 
and quadratic attraction in two dimensions. The initial data is the function 
Po(x,y ) = C(1 — x 2 — y 2 )\, with C chosen so po has mass one. Solutions are 
generated by the blob method (31) with e = 0.20,0.10, and 0.03 on the square 
[—1,1] x [— 1 , 1], which is discretized on a 80 x 80 grid. As e —> 0, the support of 
the numerical solution spreads evenly throughout the ball. This suggests that, 
in the limit, the solution converges to the indicator function on the ball, the 
unique steady state of the PDE. 


Inspired by these results, we consider a regularization of the nonlocal interaction energy 
analogous to Bertozzi and the first author’s blob method. Specifically, given a smooth, radial, 
rapidly decreasing modifier ip, we define the regularized kernel 

K e := ip e * K * ip e , 

with if e {x) := e~ d (p(x/e), and consider the regularized interaction energy 


( 1 . 6 ) 


M/*) : = 


K e (x - y) dp(x)dp(y). 


We show that, for a range of repulsive-attractive kernels, including repulsive-attractive power- 
law kernels (1.3) with 2 — d^p<0<q^2, the regularized energies E e T-converge to E with 
respect to the 2-Wasserstein metric. Then, via a compactness argument, we show that, up to 
a subsequence, the nrinimizers of E e converge to a minimizer of E as e —> 0. This provides a 
method of approximating nrinimizers of E by nrinimizers of energies with superior convexity 


and regularity properties (see Remarks 2.7 and 2.8). It also demonstrates a type of continuity 


among nrinimizers of E e , in spite of their vastly different regularity properties, as illustrated in 
Figure [lj Finally, it imparts further theoretical justification for the success of the Bertozzi and 
the first author’s numerical blob method. Our work builds on previous work of Fellner and 
Raoul 1341, who considered regularization of the Coulomb repulsion in one dimension, showing 
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that, in the presence of a confining potential, Dirac-type stationary states converge to a unique 
L°° stationary state. 

We conclude by considering the convergence of the corresponding gradient flows for repulsive- 
attractive power-law kernels (1.3) in the regime 2 — d^p<0<q^2. Indeed, using the 
abstract scheme of Serfaty 0 , we show that as e —> 0 gradient flows of E e that are bounded in 
L°°(M rf ) T-converge to a generalized notion of gradient flow for E, that is to its curve of maximal 


slope on the metric space of probability measures with bounded density. This provides a link 
between the well-understood case of gradient flows of semi-convex energies E e and the curve 
of maximal slope of the unregularized energy E. It also provides a first step in understanding 
the connection between the gradient flows of unregularized and noncovex interaction energies 


E and the aggregation equation, via a singular perturbation approach (see Remark 5.2). 

These results provide many directions for future work. A natural extension would be to 
consider the effect of simultaneously discretizing and regularizing the energy and to study the 
distinguished limits as the discretization and regularization are removed. Another interesting 
direction would be to seek quantitative estimates on the rate of convergence of minimizers and 
gradient flows as a regularization is removed. Our present work strongly leverages compactness 
arguments, and we believe that a very different approach would be needed in order to provide 
quantitative estimates. 


Remark 1.1 (Power-law kernels and generalizations). In sections [2] and [3j we choose to state 
our results for interaction kernels in the power-law form (1.3), not only because these kernels 


have attracted much interest in connection with the aggregation equation (cf. 15 32 35,36 63 


and nonlocal interaction energies (cf. |19,21,28]), but also because these kernels capture the 
crucial behavior of short-range and long-range interactions, i.e., at the origin and at infinity. 
Perhaps one of the most important properties of such kernels, which we exploit in our analysis, 
is that the attractive and repulsive parts of the kernel are given separately by relatively simple 
potentials. This allows us to quantitatively compare the attractive and repulsive parts of the 
energy. Indeed, it is due to the comparison (or rather competition) between these two terms 
that minimization of such energies exhibit energy-driven pattern formation. 

Our analysis depends mostly on the behavior of K for very large and very small values of 
|x|; hence, there are many possible ways of extending our results in sections [2] and [3] to a wider 
class of kernels. In section [4] we provide general conditions on the regularity, growth at infinity, 
and behavior around singularities of K for which our results continue to hold. In particular, 
our results hold in any dimension d ^ 1 for kernels that exhibit power-law repulsion at the 
origin and power-law attraction at infinity, including logarithmic repulsion in d ^ 2. 


Remark 1.2 (Choices of parameters in power-law kernels). In sections [2] and [3] we restrict our 
attention to repulsive-attractive power law kernels satisfying (1.3) with 2 — d^p<0<q^2 
in dimension d ^ 3. We require q ^ 2 so that the attractive term and its regularization 
are functions with at most quadratic growth and their convolution with fi G T^O^) is finitely 
valued (see Remark 2.9). We require q > 0 to obtain compactness of the sequence of minimizers 
corresponding to the regularized energies E e (see Remark 3.8). We take p < 0 since our results 
are trivially true for nonsingular attractive-repulsive power-law kernels (see section [4]). On the 
other hand, we take 2 — d ^ p in order to ensure nronotonicity of the repulsive part of the 
regularized kernels with respect to e, when the modifier is the heat kernel (see Proposition |2.6[). 
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Our paper is organized as follows. In section [2j we recall fundamental results on the 2- 
Wasserstein metric and the regularization of probability measures. We use these to establish the 


above mentioned nronotonicity property of the regularized energies ( 1 . 6 ) in the regularization 
parameter e > 0. In section [3j we prove the T-convergence of E e to E and the corresponding 
convergence of minimizers. In section [4j we extend these results to more general interaction 
kernels. Finally, in section [5j we introduce background on curves of maximal slope and prove 
the T-convergence of these generalized gradient flows for repulsive-attractive power-law kernels 
restricted to bounded densities. 

Throughout, we denote by L p (M . d ) the space of functions whose p-th power is integrable 
with respect to the Lebesgue measure on R d and by L p (p) those functions whose p-th power 
is integrable with respect to the probability measure p £ V(R. d ). C(R d ) denotes the space of 
continuous functions on M d , C 00 (M d ) denotes the space of smooth functions on M d , and B r (x) 
denotes the open ball of radius r centered at a; £ M rf . All constants which appear in the text 
are positive and may change from line to line. We express the dependence of constants to 
parameters by a subscript such as C ( j p when necessary. 

2. Preliminaries 

2.1. The 2-Wasserstein metric and regularization of probability measures. We con¬ 
sider the energies E e and E over the space 


V 2 (R d ) : = ip £ p(M d ): f |x | 2 dp(s) < +oo 

l J R d 


of probability measures with finite second moment. We endow this space with the 2-Wasserstein 
metric, which we recall briefly now, along with some of its basic properties. For further 
background, we refer the reader to the books by Ambrosio, Gigli and Savare 11] and Villani 61 . 
The 2-Wasserstein distance between p, v £ V2(R d ) is 

( 2 . 1 ) d w {n. u) := f ruin ( / [ \x - y \ 2 d^{x,y)\ ^ £ C{y,u)\] 

V Ur* Jr* )) 

where C(p, v) is the set of transport plans between p and u, 

C(p, v) := 1 7 £ V(R d x R d ): (vri ) # 7 = p and (tt 2 ) # 7 = v j . 

Here 7Ti, 112 denote the projections -K\(x,y) = x and 712 (x, y) = y. For i = 1,2, denotes 

the pushforward of 7 defined by := 7 for any measurable set U C R d . 

The minimization problem (2.1) admits a solution, i.e., there exists an optimal transport 
plan 70 £ Co(/G u) so that 


dwi^i u ) ~ 


\x~ p| 2 d7o(x,y). 


Moreover, (V 2 (R d ),dw) is a complete and separable metric space, and convergence in ( 7 ^ 2 (M rf ), dw) 
can be characterized as follows: 

diy(pn) p) —>• 0 p n —> p weak-* in P(M rf ) and f Rd |x | 2 dfi n (x) —> f Rd |x | 2 dp(x), 

f R d f{x) dn n (x) -> f Rd f(x) dfi(x), 

V/ £ C(R d ) such that |/(x)| ^ C( 1 + \x — xq| 2 ). 
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We will refer to functions satisfying |/(x)| ^ C(1 + \x — xo| 2 ), for some C > 0 and xq £ M d , as 
functions with at most quadratic growth. 

In addition to regularizing our energy functionals by convolution with a modifier, we will 
also regularize our measures. We recall the definition of the convolution of a measure. 

Definition 2.1. For p £ V 2 (R d ) and <p continuous with at most quadratic growth, p * ip is 
defined by 

(2.2) [ f(x)d(n*ip)(x) = [ f * tp(y)dp(y), 

JR d JR d 

for all bounded measurable functions f : R d —> M. 


Remark 2.2. Note that p * ip is absolutely continuous with respect to Lebesgue measure and 
d(p * ip) = p * <p{x) dx, where p * <p(x) = f Rd ip(x — y) dp{y). If, in addition, ip(x) £ 
then p * p(x) £ L°°(M rf ) and (2.2) holds for all / £ L 1 (R d ). 


We recall the following lemma on the approximation of measures by convolution. 

Lemma 2.3 (cf. 11, Lemma 7.1.10]). Fix p £ 7^2(K d ) and a mollifier ip £ C'°°(M' i ) with finite 
second moment, ip ^ 0, and f Rd ip(x)dx = 1. Then for (p e (x) := e~ d ip(x/e), p * ip e £ V 2 {R d ) 
and 

(2.3) dw(p * Tei p) ^ e ( / \x\ 2 ip(x)dx) 

\Ju d J 

2.2. Regularization of energies. In what follows, we consider mollifiers which satisfy the 
following assumptions: 

(Ml) ip £ C'°°(M d ) and ip(x) = <p{ |x|) ^ 0 for all x £ R d , 

(M2) J Rd ip(x)dx = 1 , 

(M3) ip(x) ^ C\x\~ l for some l ^ max{2d + 1,4} and C > 0. 


Remark 2.4 (Finite second moment of mollifiers). Any mollifier ip satisfying the assumptions 
(M1)-(M3) has finite second moment, i.e., f Rd \x\ 2 ip{x) dx < +oo. 

As a consequence of the previous results on regularization of measures in the Wasserstein 
metric, we have the following lemma which relates E e to E. 


Lemma 2.5. Consider an interaction energy E (p) = ff RdxRd K(x — y) dp(x)dp(y), where 
K(x) is an even, nonnegative, locally integrable function that has most quadratic growth and 
is continuous outside of some ball. Given a mollifier ip satisfying (M1)-(M3), then for all 
p£V 2 (M d ) 


E e{p) = E (p*Pe)- 
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Proof. Define the right shift operator r y : f(x) H > f(x — y ). For any / and g even, 

Ty[f * g\(x) = f * g(x -y)= f(x-y- z)g(z) dz 

JHL d 

= T y [f](x-z)g(z)dz = T y [f]*g(x). 

JS. d 

Furthermore, sinc e / * g(x — y) = / * g{y — x), we also have r y [f * g](x) = r x [f] * g(y). Since 
tp E L°° (M rf ), (2.2) holds for all continuous / of at most quadratic growth and for all / E L 1 ' 


As we may decompose K into the sum of a continuous function of at most quadratic growth 
(away from the origin) and an integrable component (near the origin), the result then follows: 


Up) = 



Pe* K * p e (x - y) dg(x)dg(y) = 
p e * r y [K] * <p € (x) dg(x)dg(y) = 
p e * t x [K] (y) d[g* tp e \ (x)dg(y) = 



T y [p e * K]* p e (x) dg(x)dg(y), 

l 

p e * t v [K](x) d[n * p e \(x)dg(y), 
K(x - y) d[g. * p t \(y)d[fi * tp € ](x), 


= HP*Pe)- 


□ 


For general interaction kernels and mollifiers, there is no uniform relation between the size 
of E e (/r) compared to E(/r). Consequently, in this section and section[3| we consider interaction 
energies of the following form: 

(El) E(/i) = ff RdxR d K(x-y) dg(x)dg(y) for K(x) = \x\ q /q - \x\ p /p, 

(E2) 2-d^p<0<q^2. 

If we mollify energies satisfying (E1)-(E2) via the heat kernel, we obtain monotonicity in e 
of the repulsive part of the interaction energy. This generalizes a property used by Blanchet, 
Carlen, and Carrillo [l7| for the Newtonian potential. In section |4j we describe how to extend 
these results to more general interaction energies. 

Let K a := | \x\ q denote the attractive part of the interaction kernel and K r := —1 \x\ p the 
repulsive part. Similarly, let E a denote the attractive part of the energy and E r denote the 
repulsive part, 

(2.4) E a (p)=f f K a (x - y)dfi(x)dp(y), E ' (g) = I f K r (x-y)dp(x)dp(y). 

JR d JlBL d JR d JS. d 


Proposition 2.6 (Monotonicity of E£). Suppose the energy E satisfies (E1)-(E2) and the 
regularization mollifier is given by the heat kernel. Then for all e\ ^ £2 > 0 and g E V r 2 (U) 


E r eM < KM- 


Proof. Since if e (x) = if(e,x) = (47re) d / 2 e l x l 2 / 4e , dpif = A x ip for all e > 0 and x E M rf . In the 
parameter regime 2 — d < p < 0, the Riesz potential 


Id+ P (f)( x ) = C Tp 


f{y)\x ~ y\ p dy = —C d)P pK r * f(x) 
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satisfies the identity 

W A /) = A I d+p (f) = -I d+p - 2 (f) 

for any Schwartz class function / (see e.g. (55]). Moreover, I d+p {f) ^ 0 for any / ^ 0. 
Defining K J (x) = ip e * K r * and using the above identity we obtain 

. K r e {x) = 2i/j e * K r * (de'ipe) = 2^ t * K r * (A'i/ig) = -(pC dtP )~Ve * I d+p (Aip e ) 

= (pC d)P )~ l 'ipt * i p+d _ 2 (A) ^ o 

since p < 0 and ip £ ^ 0. Thus, K r t is monotonically decreasing in e. 

When p = 2 — d, I d+p (f) is the Newtonian potential and satisfies the identity I d+p {Af) = 
AI d+p (f) = —/. Again, differentiating ATJ with respect to e yields the desired result. Therefore, 
by definition of E r (2.4), we conclude monotonicity in e. □ 


We close this section with brief remarks on the convexity and differentiability of regularized 
energies. 


Remark 2.7 (A e -convexity of E e ). Recall that a function / : —> M U {+oo} is A-convex if 

/(x) — ^|x | 2 is convex for some A 0. If / is twice continuously differentiable, this is equivalent 
to D 2 f(x ) ^ XI for all x E M rf . Consequently, if <p satisfies (M1)-(M3) and E satisfies (El)- 
(E2), then both the regularized potentials K e and their opposites — K e are A e -convex with 
A e = —C tp e~ d , where C^ is a positive constant depending on p and e is sufficiently small. For 


an interaction energy of the form (1.1), A-convexity of the kernel K ensures A-convexity of the 
energy E with respect to the Wasserstein metric 22 . Thus, E e and — E e are both A e -convex. 


Remark 2.8 (Differentiability of E e ). Given a functional F : V 2 {U d ) —> M U {+oo} and p E 
V 2 (M. d ) so that F(/x) < +oo, the metric local slope of F at p is 


|(9F|( p ) := limsup 


(F (p) - F(p)) + 
dw{n,v) 


If p satisfies (M1)-(M3) and E satisfies (E1)-(E2), the metric local slope for E e is well-defined 
(cf. [22} Proposition 2.2], [l] Lemma 10.1.5]) and 

\dE e \(p) = 2\\VK e *p\\ L 2 M . 

Furthermore, for p E ^(K^), we have 2||ViL e * p \\ l 2 ( p ) ^ C Vtll e 1 ~ d for e sufficiently small. 
To see this, let 0 ^ p(x) ^ 1 be a smooth function satisfying p(x) = 0 if |x| ^ 1 and 
p{x) = 1 if |x| ^ 2 . Using p, we decompose Viv e into its singular and nonsingular components, 
X7K e = X7K t ( 1 — 77 ) + X7K e p = VA' e s + V K™. By linearity of convolution and the triangle 
inequality for L 2 (p), it suffices to estimate | V Kf * p\\l 2 (ij,) and | ViL” * p\\ L 2 ( p y The former is 
bounded by ||ViL| * p\\ L°°( dp ) ^ su PxeK d |ViL|(x)| ^ C lp e 1 ~ d . To bound the latter, we apply 
Minkowski’s integral inequality, the fact that VAT™ has at most linear growth, and the fact 
that p has finite second moment to obtain 


1/2 
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Remark 2.9 (Attraction powers q > 2). One possible way to extend our results to cover 
attraction powers q > 2 would be to consider the energy E over the space V q (M. d ) of probability 
measures with finite moments up to order q, endowed with the g-Wasserstein distance 


d q (y, v) 


^min 


x - y\ q d^(x,y): 7 <E C{y,u) 


1/9 


where C(y, v) is the set of transport plans between y and v. Alternatively, to cover all q > 0, 
one could consider the oo-Wasserstein metric 


doo(y, u) := inf sup \x - y\ 

7(i,j/)£supp( 7 ) 

over the space Voo(M d ) of probability measures with finite moments of all orders. The space 
{V oc ), doQ ) is a complete metric space 


37 


and has been used in the study of local minimizers 
for several variational problems, including nonlocal repulsive-attractive energies |6,21,46 . 

The main difficulty in extending our results to V q (M. d ) or V oa (W i (ii) (iii) ) lies in understanding 
the appropriate generalization of the A-convexity to these distances and, in particular, the 
appropriate analogue of the HWI inequality, which plays a key role in the proof of Theorem 
13.41 We leave this to future work. 


3. Convergence of regularized energies and minimizers 

We now turn to the proof that, up to a subsequence, minimizers of the regularized energies 
E e converge to a minimizer y of E with respect to d\y- We establish this by first proving a 
T-convergence result for the sequence of energies E f and then obtaining compactness for any 
sequence {y e } e> 0 when E e (/r e ) is uniformly bounded. For the latter, we use Lions’ result on 
concentration compactness, which we recall for the readers’ convenience. 


Lemma 3.1 (Concentration compactness lemma for measures (cf. [43], [56| Section 4.3])). 
Let {y n }neN be a sequence of probability measures on M rf . Then there exists a subsequence 
{Rn k }keN satisfying one of the three following possibilities: 


(i) (tightness up to translation) There exists a sequence {ykfkeN C 
e > 0 there exists R > 0 with the property that 


such that for all 


' B R (yk) 


dy nk (x) ^ 1 — e for all k. 


(ii) (vanishing) lim sup / dy nk (x) = 0, for all R > 0; 

fc_s>0 ° y<m. d JB R (y) 

(iii) (dichotomy) There exists a E (0,1) such that for all e > 0, there exist a number R > 0 
and a sequence {x/JfcgN C with the following property: 

Given any R' > R there are nonnegative measures y\ and yj, such that 

0 < Mfc + Mfc < hn k , supp(/4) c B r ( x k ), supp (yl) cm d \ B RI {x k ), 


lim sup 

k —^00 


dyl{x) 


+ 


(1 - a) 



^ e. 
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We begin our proof of the convergence of minimizers with the observation that minimizers 
exist for both the regularized and unregularized energies. 


Proposition 3.2 (Existence of minimizers in T^Q^))- Suppose the energy E satisfies (El)- 
(E2) and the mollifier ip satisfies (M1)-(M3). Then both E and E e admit a minimizer in 
7^2 (M rf ), for e > 0 sufficiently small. 


Proof. First note that, in the parameter regime 2 — d<p<0<q^2, the interaction potential 
K and its regularization K e are locally integrable and lower semicontinuous. 

We now show that K and K e are strictly increasing in each coordinate for Xi sufficiently 
large. For K given as in (E1)-(E2) and x 0, 

d Xi K(x ) = Xi\x \ q ~ 2 — Xi\x \ p ~ 2 = Xi\x\ q ~ 2 (l — \x\ p ~ q ). 

Consequently, if x% > 1, d Xi K(x ) is positive. 

Now we consider K e . It suffices to show that if Xi > 2, then 

(3.1) \d Xi K e (x) - d Xi K{x)\ < Ce\x\ q ~ 2 


for some constant C > 0. Then for Xi > 2, we have (1 — \x\ p q ) > C p . q > 0 and 

d Xi K e (x) = d Xi K(x ) + d Xi K e (x) - d Xi K(x) Xi\x\ q ~ 2 {l - |ic| p_ ' ? ) - Ce\x\ q ~ 2 

> 2\x\ q ~ 2 C Py q - Ce\x\ q ~ 2 = \x\ q ~ 2 {2C P) q - Ce). 

Choosing e sufficiently small, this is nonnegative. 

We now prove inequality (3.1). First, we use a cutoff function to rewrite d Xi K(x) as the sum 
of a compactly supported singular function d Xi K s and a continuously differentiable function 
d Xi K n . Let 0 ^ i/ ^ 1 be a smooth function satisfying rj(x) = 0 for |x| < 1/4 and r](x) = 1 
for |x| > 1/2. Write d Xi K = (1 - p)d Xi K + r]d Xi K =: d Xi K s + d Xi K n . Then, for <3? e (y) := 
p e (y- z)ip e (z)dz, 


| d Xi K e {x) - d Xi K{x)\ = 


{■d Xl K(x - y) - d Xi K{x)) <& e {y) dy 




[ (d Xi K s (x-y)-d Xi K s (x))$ e (y) 
JR d 


dy 


+ 


(d Xz K n (x - y) - 8 Xi K n {x )) 4> e (y) dy 


=:h+h- 


The autocorrelation function 4> e satisfies 4 > e (i/) = e~ d &(y/e) for all y € and has the same 
decay property (M3) as ip e . To see this, suppose |y| = 3 R for some R > 0. Then for any z G 
we have \z\ > R or \y — z\ > R] hence, 



< Ce l ~ 2d R~ l = Ce l ~ 2d \y\~ l . 
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II 


First, we estimate I If Xi > 2, then |x| > 2, so d Xi K s (x) = 0. Consequently, 

h = [ d Xi K s (x-y)$ e (y)dy . 

J\x-y\<l 

Since \x — y\ < 1 and |z| > 2, 

\y\ > M - \y - x\ > |x| - 1 ^ \x\/2. 

Consequently, 

|^(y)| ^ Ce l ~ 2d \y\~ l < Ce\y\~ l < C2 l e\x\~ l 

for 0 < e < 1. Since d Xi K s is locally integrable and 2 — q < 2d + 1 ^ l, there exists C > 0 so 
that for |x| > 2, 

h ^ Ce \x\- 1 ^ Ce|x| 9 - 2 . 

Now we estimate / 2 - Since d Xi K n is continuously differentiable, 


h = 



1 d 

—d Xi K n {x - ay) da $ e (y) dy 

q CiO!. 

1 

{Vd Xz K n (x - ay), -y) da $ e (y) dy 




R d Jo 


1<l*l/2-/0 


\Vd Xi K n (x - ay)\ \y\da $ e (y)dy 


u 


+ / / \Vd Xi K n (x - ay)\\y\da$ e (y)dy. 

J \ y \>\ x \/ 2 jo 

For \y\ ^ \x\/2, \x — ay\ ^ \x\/2 > 1 for all a 6 [0,1]. Thus, V7y(x — ay) = 0, and 
\S7d Xi K n {x - ay )| ^ | (Vrj{x - ay))d Xi K(x - ay )| + \rj(x - ay)Vd Xi K(x - ay)\ ^ C\x\ q ~ 2 . 
On the other hand, for all all y G M d , |V9 a , i lF ri (x — ay) | ^ C. Therefore, 

h^C\x \ q ~ 2 ( \y\<S> e (y)dy+ C [ \y\$> € (y)dy 

J\v\<\x\/2 


l>M/2 


= Ce|x| 9 2 f \z\ <h(z) dz T Ce [ \z\®(z)dz. 

j\z\<\x\/(2e) J\z\>\x\/{2e) 

The first integral is bounded by a constant since <1? satishes decay property (M3), hence, has 
finite first moment. For the second integral, we use that |z<l>(, 2 ;)| ^ C\z\~ l+1 and 2 — q< l — d— 1 
to obtain 


<\z\>\x\/(2e) 


|z| ^(z) dz 


r+oo 


r -i+l r d-l dr = C ( p. ] ^ Ce f x 


l\x\/(2e) 


2e 


d-l+l 


I < 7-2 


Therefore, (3.1) follows and we conclude that both K(x) and K f {x) are strictly increasing in 
each coordinate for |x| > 2. 

We now show that K{x) and K e (x) become arbitrarily large as |x| -A Too. This is immediate 
for K(x), since q > 0. Furthermore, since K(x) and $ e (ic) are nonnegative and f Rd <S> e (x) dx = 
1, we also have 

|x|—>-+oo 

i\\x — y)^ E \y)ay ^ mm i\{Z) 

KN /2 


K e (x) ^ / K(x — y)<& e (y)dy ^ min K(z) 

.1 V x 12 N5»M/2 


A Too. 










12 


KATY CRAIG AND IHSAN TOPALOGLU 


Thus, for both E and E e , there exists a minimizer in [54[ Theorem 3.1]. This minimizer 

is compactly supported when e > 0 is sufficiently small, thus it belongs to 7^2 (J£ d ) [l8, Theorem 
1.4]. □ 

Remark 3.3. Since C V(R d ), the infimum of the energy E or E e over R(IR rf ) is less 

than or equal to the infimum over V 2 (K d ). Due to the fact that minimizers have compact 
support |18[ Lemma 2.10], the converse holds, as well. Hence, 

inf E(/i) = inf E(/t) and inf E e (//) = inf E e (/r) 

fi£V(R d ) 2(K d ) fi£V(R d ) iieV2m d ) 

for all e > 0 sufficiently small. 


Now we prove the regularized energies E e converge to E in the sense of T-convergence. 


Theorem 3.4 (T-convergence of regularized energies). Suppose the energy E satisfies (El)- 
(E2) and the mollifienp satisfies (M1)-(M3). Then the sequence of regularized energies { E e } e> o 
r -converges to the energy with respect to (V 2 (R rf ),dyv)- That is, 

(i) (Lower semicontinuity) For any {p e } e >o C ^(U^) and p £ R 2 (R d ) such that it holds 
lim e _^o dwihe, p) = 0, we have that 

liminf EJp e ) E (p). 

e —>0 

(ii) (Recovery sequence) For any p £ 7^2 (® d ) there exists {Pe}e >0 c 7*2 Q^) such that 

lim dw(v e , p) = 0 and lim E e (i/ e ) = E (p). 

e ->0 e —>0 


Proof. We will prove this theorem in two steps. 


Step 1 (Lower semicontinuity). By Lemma 2.5, we have E e (p e 


E (p e *p e ). Furthermore, 


(3.2) 


d W {p, Pe * Ft) ^ d W (p, Pe) + dw(Pe, he * Te), 


where the first term approaches zero by hypothesis and the second term approaches zero by 
Lemma 12.31 

As the interaction potential K is lower semicontinuous and bounded from below, the Port¬ 
manteau Theorem |60[ Theorem 1.3.4] ensures the energy E is lower semicontinuous with 
respect to convergence in V 2 (M d ). Indeed, the Portmanteau Theorem states that the weak-* 
convergence of p e x p e to p x p is equivalent to the fact that 


lim inf 

e —>-0 



y)d(p e X Pe){x,y) 



y)d{p X p)(x,y) 


when / is lower-semicontinuous and bounded from below. Hence, we obtain 


liminf EJp e ) = liminf E(p e * ipf) E (p). 

e ->0 e ->0 


Step 2 (Recovery sequence). Let p £ T^O^) be arbitrary. We define a recovery sequence 
for the measure p by using the heat kernel. Let ip(x) = (47r) -d//2 e~^ / 4 , and define 
fis(e)(x) '■= ^( e )~ d '4’(x/6(e)) where 


(3.3) 


1/2 d 


6(e) := e 
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Clearly the function if) satisfies the assumptions (M1)-(M3) in Subsection 2.2 hence, it is an 
admissible modifier. Define the measure 


Then, by Lemma 


2.3 


v e := M* V’S(e)- 

v e G T , 2 (K d ) for all e > 0, and dw{v e , r) -A 0 as e 


0 . 


We now show that lim e _>o E e (zz e ) = E(yn). By Lemma 2.5 


(3.4) 


|E e (u e ) - E(/x)| = |E(/X * V’S(e) * ¥>e) - E(/i)| = |E 5(e) (/i * <p e ) - E(/x)| 


< I E <5(e) (/^ * ¥>e) - E (5(e ) (//,)| + |E 5(e) (/i) - E(/i)| =: h + I 2 . 

It suffices to show that I \, L> —> 0 as e —> 0. We estimate I 2 first. As E ste)^) is regularized 
using the heat kernel ips(e) > Proposition 


2.6 


ensures 


that E^(aO ^ E r (^) for all e > 0. Hence, 


limsup E 5 ( e )(/z) ^ E r (/i) + lim sup E£ (e) (/i) 


e —>0 


e-s>0 


Again, by Lemma 2.5, E^(//) = E a (/r * Since dw(R * i/ts( e ), /i) —> 0 as e —> 0 and the 


attractive interaction kernel K a is a continuous function with at most quadratic growth, 


Therefore 


limsup E% e) (n) = limsup E a (/r * V’ff(e)) = E “ (/■»)■ 

e —>0 e —>0 


limsupE 5(e) (^) ^ E(/i). 

e —>0 


On the other hand, by the lower semicontinuity of the energy E, as proved in the first step of 
the proof, we get that 

iinDnf E S (e)(v) = liminf E(/x * ^(e)) ^ E (/h>; 
hence, I 2 = |E,j( e )(/r) — E(pi)| -A 0 as e -A 0. 

To estimate 1 1 , recall that by Remark |2.7| and [2.8| the regularized energies ±E d -( e ) are \s(e)~ 
convex with X s ^ = —C lf 5(e)~ d and have metric local slope ||VA * i!>s{e) * pIIl 2 ^)- Therefore, 
they satisfy the following HWI-type inequality jl, Theorem 2.4.9], 

I E <5(e)(M*¥k) - E 5(6)(/L)| < 2|| VH*4( e ) * ^\\ L 2 M d W (^* ^Pe,^) -* ip e , (J,). 


Consequently, for e > 0 sufficiently small, so that dw (fj, * ip e , //) ^ 1, we may use the bound 
||VK * * A 4 11 L 2 (fj.) ^ C^^5{e) 1 ~ d f rom Remark 2.8 and the definition of 5(e) from equation 


to obtain 

h = I E 5 (e) (/i * <pe) - E tf(e) (/i)| < ^(e) 1 ^ + 5(e)~ d ) d w (n * <p e ,n) < Ce5(e)~ d = Ce 1 ' 2 
Therefore I\ -A 0 as e —> 0, which completes the proof. 


□ 


Remark 3.5. Since the energy E is lower semicontinuous with respect to weak-* convergence 
of measures in V(M. d ) and convergence in V 2 (M. d ) implies weak-* convergence in V(M. d ), the 
conclusion of Theorem 3.4 4) is also true with respect to weak-* convergence in T’(M d ). 


As a result of the strong confining forces induced by the attractive part of the kernel K, we 
obtain the following compactness result for the sequence of energies E e in P(M rf ). 
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Proposition 3.6 (Compactness in P(M rf )). Suppose the energy E satisfies (E1)-(E2) and the 
mollifier ip satisfies (Ml)-(M3). Let {y e } e> o C V(M. d ) he a sequence such that, for all e > 0 
sufficiently small, E e (y e ) ^ C for some constant C > 0. Then {y e } e> o has a subsequence which 
is, up to translations, convergent with respect to the weak-* topology in V(R d ). 


Proof. We will use Lemma 3.1 and argue by contradiction, as in 54 . Note that as in the proof 
of Proposition 3.2, K(x) and K e (x) —> +oo as |x| —> +oo for e > 0 sufficiently small. This 
ensures that the energy E e exhibits long-range confinement. We will use this fact to eliminate 
“vanishing” and “dichotomy” possibilities of the sequence {y e } e> o. 

Suppose a subsequence of {y e } e> o, which we still index by e > 0, “vanishes” in the sense of 
Lemma |3.1[ ii). Then for any 6 > 0 and R > 0, there exists eo > 0 such that for all e < eo and 
i 6 K d , we have 

li e (R d \ B r (x)) >1-5. 

Hence, for all e < eo, 


(3.5) / / dy e (x)dy e (y) = / / dyfiy) dyfix) >1-5. 

JJ\x-y\^R JR d \JR d \B R (x) ) 

Since K e (x) —> +oo as |x| —> oo, given a constant B > 0, there exists R > 0 such that for all 
|x| > R, we have K e (x) > B. Let 5 < 1/2, and choose eo depending on 5 and R such that 
(3.5) holds. Then, using (3.5) and the fact that K e > 0 for all e > 0, we get that 


= 0-0 = // _K e (x-y)dy t (x)dy e (y) + _ K e (x - y) dy e (x)dy e (y) 

J J\x—y\<R J J\x—y\^R 


> 


\x-y\ <R J J\x-y\^R 

K e (x - y) dy e (x)dy e (y) > (1 - 5)B. 

\x-y\^R 


This contradicts the uniform bound E e (y e ) Sj C when B is sufficiently large; hence, “vanishing” 
does not occur. 

Suppose that for a subsequence of {y e } e> o, which again we index by e > 0, “dichotomy” 
occurs. Then, since K e > 0 for all e > 0, 


liminf E e ( / u e ) > liminf / / 

e ^° e ^° JR d \B' R {x,)J B R {x e ) 

> inf KJx)a( 1 — a), 

| x\2R'-R 


K e (x - y) dy\(x)dy 2 e (y) 


where R > o, the sequence {x e } e >o, and the measures y\ and y 2 are defined as in Lemma 
|3.1[ iii), and R' > R is arbitrary. Thus, again using the growth of K e , we get that 

liminf EJy e ) > +oo, 

e ->0 

contradicting the uniform bound on E e (^ e ). 

Therefore, Lemma 3.1 ensures the sequence {y e }e>o is tight up to a translation by a sequence 
{y e } e >o C R d . However, the energies E e are translation invariant, so we can replace y e by 
y e f + y t ) which yields a tight sequence. Therefore by Prokhorov’s theorem (cf. [60, Theorem 
1.3.9]), there exists a subsequence {y e } e> o, still indexed by e > 0, and a measure y S V(R d ) 
such that, up to translations, y e y with respect to weak-* convergence in V(R d ) as e —> 0. □ 
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Classically, if a sequence of T-convergent functionals also satisfies a compactness property, 
then, up to a subsequence, minimizers co nverg e to a minimizer of the limiting functional. 
Though our compactness result, Proposition 3.6, is established in a weaker space (V(M. d )) than 
the topology in which the energies T-converge (TL^R^)); we still obtain the following corollary 
on the convergence of minimizers. 

Corollary 3.7 (Convergence of minimizers). Suppose the energy E satisfies (E1)-(E2) and 
the mollifier satisfies (M1)-(M3). For any e > 0 sufficiently small let g £ £ V 2 (M d ) be a 
minimizer of the energy E e . Then there exists /a £ such that, up to a subsequence and 

translations, g £ —> g in as e —> 0, and g minimizes the energy E over T^O^)- 

Proof. Since the sequence {g £ } £> o C P 2 (R d ) consists of minimizers of E e , for e > 0 sufficiently 
small, there exists a constant C > 0 so that Efigf) ^ C. Then by Proposition 3.6, there exists 
a measure g £ V(M. d ) such that a subsequence of {g £ } £> o, which we still denote by g £ , weak-* 
converges to g in V(M. d ). 

To show that g minimizes E over we proceed in two steps. First, consider an arbitrary 

measure v £ By Theorem 3.4 (ii), there exists a sequence {u e } e >o C so that 

dwffie-, u) —^ 0 as e —y 0 and 

limE(u e ) = E(u). 

e-s>0 


By Theorem 3.4 (i), Remark 3.5, and the fact that the measures g £ are minimizers of E e over 
V 2 (R d ) C 'P(M 3 ), we obtain 

(3.6) 

Thus, g £ 


E (g) ^ liminf EJg e ) ^ liminf E e (zz e ) = E(u) 

e —>0 e — 

( ) has less energy than any other measure u £ so by Remark 


3.3 


18, Lemma 2.10] the authors 


g minimizes E over V(M. d ), as well. Consequently, |l8j Lemma 2.10] ensures g is compactly 
supported; hence it is in 

Finally, we show that in fact g £ converges to g in T^O^)- In 
show that 

diam(supp g e ) ^ M £ 

where M e := \/d(Ar £ + (|"1 /m £ ~\ — l)(4r e + 2 Rff). Here 

C-E £ (g £ ) 


m f := 


C - 


with C > 0 so that E e (g t ) ^ C for e > 0 sufficiently small. Kfi 110 denotes the absolute minimum 
value of the interaction kernel K £ , r e is chosen such that K e (x) C for all |x| ^ r e , and R e 
is the radius after which the kernel K £ is strictly increasing. We denote the corresponding 
quantities for the unregularized energy E by removing the subscript e. 

We show that M £ , the upper bound on the diameter of the support of g e , is bounded by 
some M > 0 for all e > 0 sufficiently small. We have limsup e ^. 0 — E £ (g £ ) ^ — E (g) and 
K e — > K uniformly on compact sets away from the origin. Thus, m £ — > m as e —> 0, so 
\l/mf\ ^ |T/m] + 1 for e > 0 sufficiently small. 

As shown in the proof of Proposition 3.2 K e is strictly increasing for |ir| ^ 2, so we may 
take R e = 2. Likewise, by estimate ( |3.1 ), there exists r > 2 so that K e (x) C for |x| > r. 
Therefore 

M e ^ M := Vd(4r + (|"l/m])(4r + 4)) 
for all e > 0 sufficiently small. 
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This shows that the diameter of the support of // e is uniformly bounded by M. Consequently, 
the second moments of the sequence y, e are uniformly integrable 


lim / 

k^+oo J{\ x \2^k} 


\x\- dy, t = lim 


fc — S >+00 




\x\ 2 dpL e = 0 


Since p e converges to p with respect to the weak-* topology on V(M. d ) and [A e has uniformly 


integrable second moments, —> /i in 7^*2(IR rf ) with respect to d\y 111, Proposition 7.1.5]. □ 


'(M 

0 


Remark 3.8 (Negative attraction power). The condition q > 0 is a sufficient condition for 
Propositions |3.2| and 3.6 If q < 0, both K and K e converge to zero as |x| -A oo. In 54 


the authors characterized the existence of minimizers for these types of potentials using the 
existence of a measure for which the energy is negative. This is clearly true for E and E e , as 
one can consider the characteristic function of a sufficiently large ball Hr( 0). However, we 


still require q > 0 for the compactness result in Proposition 3.6, and we believe that it is a 
necessary condition, as well. 

4. More general interaction energies 

Our results of sections [2] and [3] can be extended to more general nonlocal interaction energies 
(4.1) E (n)= f [ K(x-y)d/j,(x)dn(y) 


that are not necessarily repulsive-attractive and are defined via an interaction kernel K : M d — > 
MU {Too} that satisfies the following hypotheses for some fixed R > 1: 

(HI) K is even, i.e., K{x) = K(—x) for all x £ M d . 

(H2) K £ Ll oc (R d ) n C ,1 (M d \ B r ( 0)) and |VK(x)| ^ C(1 + |x|) for |s| > R. 

(H3) For |x| > R, K is strictly increasing in each coordinate and 

lim K(x) = Too. 

|x|—S>+00 

(H4) There exists a function K a £ C(M d ) so that K r := K — K a is super harmonic. 

(H5) For |x| large, K r and K a have at most quadratic growth, i.e., 

\K a (x)\ T \K r (x)\ ^ C( 1 T |x| 2 ) for |x| > R. 

We briefly recall the definition of superharmonicity and some immediate properties of su- 
perharmonic functions. 


Definition 4.1. A function f £ 

/(*) > 

for a.e. x £ 


is superharmonic if 
f(y ) dy 


B r(x)\ J Br (x)' 

and every r > 0, where \B r {x)\ denotes the volume of the ball. 


Remark 4.2 (Properties of superharmonic functions). If / £ L 1 1 oc (M‘ i ) is superharmonic, then / 
is lower semicontinuous, it is bounded below on compact sets, and its distributional Laplacian 
satisfies A/ ^ 0 [42| Theorem 9.3]. 
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Interaction energies (4.1) defined via kernels satisfying assumptions (H1)-(H5) cover a wide 
range of applications. In particular, in any dimension d ^ 1, they include energies with K of 
one of the following forms: 

• Attractive-repulsive power-law kernels ( |1.3[ ), for 2 — d^p<0<q^2. (This includes 
Newtonian repulsion, p = 2 — d, in d ^ 3.) 


Nonsingular attractive-repulsive power-law kernels (1.3), for 0 < p < q ^ 2. (This 
includes Newtonian repulsion, p = 1, in d = 1.) 

For d ^ 2, attractive-repulsive kernels with logarithmic repulsion, K(x ) = \x\ q /q — 
log(|a;|) for 0 < q ^ 2. (This includes Newtonian repulsion in d = 2.) 

Attractive-repulsive kernels with power-law behavior at the origin and infinity, i.e., 
radial kernels for which there exist 0 < R\ < R 2 and C\, C 2 > 0 so that for either 
2 — d^p<0<q^2or0<p<q^2, 

K{x) = Ci \x\ p if d ^ 1 or K(x) = — C\ log(|x|) if d ^ 2 on 13^(0); 


K{x) is continuous on Br 2 (0) \ Br a (0); and, 
K(x) = C 2 \x\ q on R d \ Br 2 ( 0). 


Remark 4.3. Any function K satisfying the hypotheses (H1)-(H5) is locally integrable, bounded 
from below, and lower semicontinuous. Consequently, without loss of generality, we will assume 
that K is nonnegative, since the minimizers of an interaction energy E with kernel K are the 
same as the minimizers of E with kernel K — inf K. 


Our results continue to hold for K satisfying (H1)-(H5), provided that we define the regu¬ 
larized energies E e via compactly supported mollifiers. That is, we replace mollifier assumption 
(M3) by 

(M3’) ip is compactly supported in M. d . 


This condition eliminates possible oscillations in K e for large values of |x|, allowing us to 
establish growth and monotonicity properties of the regularized kernels when |x| is sufficiently 
large. 

We now describe how our previous results naturally extend to K satisfying (H1)-(H5). 
The key result in section [ 2 ] is Proposition 2.6 Assumption (H4) ensures that the potentially 
singular part of the kernel, K r , is locally integrable and superharmonic. Hence, Proposition 


2.6 continues to hold since, again taking 'i/y to be the heat kernel, 


^-K r e = 2*p e * K r * (d e i/> e ) = 2ifj e *K r * (Ay> e ) = 2ift e * (A K r ) * ^ 0, 

by the superharmonicity of K r . 

The remaining results of section [ 2 ] also continue to hold for energies satisfying (H1)-(H5) 
and mollifier s sat isfying (Ml), (M2), and (M3’). In particular, (H2) ensures K e L ] 1 oc (M fi ) so 
that Remark 


2.7 


holds with A e = C^e d 2 


, since for 0 < e < 1, 


\D 2 K f 


L°° 




I D 2 p e 


\L°° 


K{x) dx ^ C v e 


—d—2 


(4.2) 


supp Ip 


< + 00 . 
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Likewise, assumption (H2) ensures that Remark 2.8 holds. As before, we may choose a cut¬ 
off function p to separate K into a locally integrable component K s near the origin and a 
continuously differentiable component K n away from the origin. However, since for 0 < e < 1, 


(4.3) 




L°° 


< ||V^ e 


\L°° 


J 

J su 


K s (x ) dx ^ C p 


— l—d 


< supp ip 

we merely obtain the bound 2||VR e * pWl 2 ^) ^ C v ^e~ l ~ d on the metric local slope. 

Now, we describe how to extend the results from section [3j Proposition 3.2 merely requires 
that K and K e are locally integrable, lower semicontinuous, strictly increasing in each coordi¬ 
nate outside of some ball, and become arbitrarily large as |x| -A +oo. To see that the last two 
properties hold for K e , note that for 0 < e < 1, 

K e (x) - K e {y ) = I [K(x - z) - K(y - z)\y e (z) dz , and 

J supp ip 
K e (x) > 


min K(x — z). 

z£ supply 


Consequently the properties for K e are immediate consequences of the corresponding properties 
for K, which hold by (H3). 

The proof of Theorem |3.4| is nearly identical for our more general energies, though we must 
take <5(e) = e 1 / 2 ( rf + 2 ) to compensate for the inferior bounds on the convexity constant A e (4.2) 
and the metric slope (4.3). Finally, the proofs of Proposition 3.6 and Corollary 3.7 again use 
that K and K e become arbitrarily large (uniformly in e) as |x| —> Too, are strictly increasing 
in each coordinate outside of some ball (uniformly in e), and are both nonnegative. 

Consequently, though we choose to state our main results in sections [2] and [3] for energy 
functionals with repulsive-attractive power-law potentials satisfying (E1)-(E2) (see Remarks 
0 and |1.2[ ) , these results naturally extend to a wider class of energy functionals 0 with K 
satisfying (H1)-(H5). 

Remark 4.4 (More singular kernels). The main difficulty in extending our results to interaction 
kernels that are more singular than the Newtonian potential at the origin (i.e., —d < p < 2 — d) 
lies in the construction of a recovery sequence in the proof of Theorem 3.4 Although one could 


obtain a monotonicity result as in Proposition |2.6| for such kernels by regularizing them via 
fractional heat kernels, unfortunately, to our knowledge, in the parameter regime —d < p < 
2 — d, the decay estimates on fractional heat kernels do not guarantee that they have finite 
second moments. Therefore they are not admissible candidates in constructing our recovery 
sequences. 


5. Convergence of gradient flows 


In this section, we apply our result on the T-convergence of regularized energies satisfying 
(E1)-(E2) to show that if the Wasserstein gradient flows of E e are bounded in L°°(M rf ), they 
converge to a metric space generalization of the gradient flow of E, known as a curve of maximal 
slope. 

As previously mentioned, without regularization, energies E satisfying (E1)-(E2) are not 
convex (or A-convex for A < 0). Consequently, they fall outside the scope of much of the 
existing theory on well-posedness of Wasserstein gradient flows [I 22]. However, the regularized 
interaction energies E e are A e = C lfi e~ d convex, so that for fixed e > 0, their Wasserstein 
gradient flows exist and are unique [22 . Thus, by showing these gradient flows T-converge 
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to the curve of maximal slope for the unregularized energy, we provide a link between the 
well-understood case of convex gradient flow and emerging results on the Wasserstein gradient 
flow of non-convex energies [2j[l5|[23|[24|[26|[30j|44j|45] . 

We restrict our attention to the space of probability measures with bounded density 


(5.1) 


V 2 , R (R d ) :={M6M Rd ) : 


L°° 




for any R > 0, due to the fact that, though E is not convex (or A-convex for A < 0), once it 
is restricted to V 2 .n(^ d )- it possesses a generalized notion of convexity known as ui-convexity, 
where uj{x) is a log-Lipschitz modulus of convexity. 


Definition 5.1 (cf. [23[|30|). Consider the modulus of convexity 

{ xllogd if 0 ^ x ^ e~ l ~^, 

r-^---— 

y / x 2 + 2(l + ^/2)e- 1 -Ar if x > e 1 

Then E : 'P 2 ,R.(M d ) —> M U {+oo} is w-convex if, for all constant speed geodesics p, a : [0,1] - 
'P'ljifpd), there exists C > 0 so that 

(j 

E (no) < (1 - a)E(/i 0 ) + aE(/ii) + — [(1 - a)u;(a: 2 d^(/zo, /h)) + au(( 1 - a) 2 d^(/i 0 , /xi))] • 


This notion has been used, explicitly and implicitly, in many previous works on non-convex 
gradient flow 23,24. For our purposes, cu-convexity allows us to define a notion of upper 


gradient for E, which is an essential component in defining its curve of maximal slope. 

Both from the perspective of energy minimization and the dynamics of the aggregation 
equation (1.4), the restriction to bounded densities is quite natural when E is a repulsive- 


attractive interaction energy satisfying (E1)-(E2) with p = 2 — d. Indeed, several results in the 
literature support the fact that both the energy minimization and evolution take place in the 
space of bounded functions. For quadratic attraction (i.e., for q = 2), the global minimizer of 
E is the characteristic function on a ball [21,28 . Likewise, from the dynamics point of view, if 


the initial data belongs to 'P 2 {^ d )C\L°°(R. d ) and has compact support, then it remains bounded 


for all time and converges to the characteristic function on the unit ball 14 . For more general 
attraction powers (0 < q ^ 2), all compactly supported local minimizers are bounded 21 and 


smooth, compactly supported classical solutions of the aggregation equation remain bounded 
for all time (5, Lemma 1]. 

For p / 2 - d, the restriction to measures with bounded density is perhaps less natural, 
since the minimizers may concentrate mass on sets of measure zero and classical solutions to 
the aggregation equation can approach these steady states asymptotically [4j. However, since 
our proof regarding the T-convergence of gradient flows works for any choice of parameters 
2 — d^p<0<q^2, we choose to include this regime for the sake completeness. 


Remark 5.2 (Gradient flow vs. aggregation equation). While the aggregation equation (1.4) 


does inform our choice of the space of measures with bounded density (5.1), the relationship 


between weak solutions of the aggregation equation and the Wasserstein gradient flow of the 
interaction energy is purely formal for the unregularized (hence nonconvex) interaction energy. 
There is hope this relationship can be made rigorous following the approach of Ambrosio and 
Serfaty [2], but we leave this for future work. 
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For the regularized energy, if g{t) £ is a weak solution of the aggregation equa¬ 

tion ( |1.4[ ), then it is also a Wasserstein gradient flow of E e [lj Corollary 11.1.8]. Below, we 
consider the curve of maximal slope of E e on the metric space of measures with bounded 
density (V 2 ,R$t d ), dw), and if g(t) £ T^O^) is a weak solution of the aggregation equa¬ 
tion with ||/i(i)||oo ^ R, then it is a curve of maximal slope for the weak upper gradient 
g e {g) := 2 \\VK e * n\\ L 2^ on (V2.n(^ d )- d\y) JlJ Theorem 11.1.3]. The reverse implication is in 
general false, due to the height constraint imposed by V 2 ,R$& d )- 


5.1. Curves of maximal slope and T-convergence. We now briefly recall the notion of 
curves of maximal slope on a compete metric space ( S,d ). We refer the reader to the book 
by Ambrosio, Gigli, and Savare 1, Chapter 1] for further details. A curve u(t) : (a,b) —> S is 
2-absolutely continuous if there exists m £ L 2 (a,b ) so that 

(5.3) d(u(t),u(s)) ^ J m(r)dr for all a < s ^ t < b. 

For any 2-absolutely continuous curve, the limit 

d(u(s),u(t )) 


|u'(t)| = lim 


s—yt 


I s — 1| 


exists for a.e. t £ (a, b ). F urthermore m{t) := |?/(t)| £ L 2 (a,b ) satisfies (5.3) and for any 
m £ L 2 (a,b ) satisfying (5.3), we have 

|i/(f)| ^ m(t) for a.e. t £ (a, b). 

Given a functional F : S -A (—oo, +oo] that is proper , i.e., D{ F) = {« £ S : F(it) < +oo} / 0, 
its upper gradient is a generalization of the modulus of the gradient from Euclidean space. 
Specifically, g : S —> [0, +oo] is a strong upper gradient for F if for every 2-absolutely continuous 
curve u(t) : (a, b) -A S the function g o u is measurable and 

(5.4) |F(u(i) — F(u(s))| ^ J g(u(r))\u'\(r) dr for all a < s ^ t < b. 

One example of a strong upper gradient is given by the metric local slope 

(F(«) - F(u))+ 


|<9F|(«) := lim sup 


d(u, v ) 


when F is a A-convex and lower semicontinuous functional 1, Corollary 2.4.10]. 

Finally, we recall the definition of a curve a maximal slope. A locally 2-absolutely continuous 
curve u : (a, b) —> S is a curve of maximal slope for F with respect to the strong upper gradient 
g if there exists a non-increasing function cj) so that <f>(t) = F o u(t) for a.e. t £ (a, b) and 


(5.5) 


< ~\\ u '\ 2 W “ \s 2 {u(t)) for a.e. t £ (a, b). 


For all e > 0, [I, Corollary 2.4.12] ensures that if g £ D( E e ), then there exists a curve of 
maximal slope g e (t): (0, +oo) -A T^rO^) for E e with respect to the strong upper gradient 
g e {y) = 2|| VA £ * z/|| L 2 ( dl/ ) satisfying g e (0) = g. 

With these definitions in hand, we now recall a general result of Serfaty on the T-convergence 
of gradient flows on a metric space. 
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Theorem 5.3 (cf. |53[ Theorem 2]). Let F e and F be functionals defined on metric spaces 
(. S e ,d e ) and ( S,d ) with strong upper gradients g e and g, respectively. Suppose the following 
criteria hold: 

(i) (T-liminf convergence) There is a notion of convergence S of u e G S e to u G S so that 

u e u implies liminf F e (u e ) ^ F(it). 

e ->0 

s 

(ii) (Lower bound on metric derivatives) Ifu e (t) —“■ u(t) fort G (0,T), then for s G [0, T), 

/*S /*S 

liminf / 

£ ->° Jo Jo 

S' 

(iii) (Lower bound on slopes) If u e —“■ it, f/ien lirn inf t _,o Se( u e) ^ g(u). 

If u e {t) is a curve of maximal slope on (0, T) /or F e with respect to g e satisfying 

uAt ) ii(f) /or f G (0, T) and lim F e (it e (0)) = F(it(0)), 

£->0 

then u(t ) is a curve of maximal slope for F with respect to g and 

lim F Ju e (t)) = F (u(t)) for all t G [0, T), 

e-S’O 

g e (u e ) -A- 5 (it) and |it'| d(E -A \u'\ d in L( oc (0,T). 

Remark 5.4. Although the metric d induces a natural topology on S, the above result admits 

5 

a notion of convergence S that can be induced by a weaker topology c on 5. Indeed, u e —^ it 
means that 7r e (u e ) u for some map 7r e : S € -A S (see 53 and |l, Remark 2.0.5] for details). 


5.2. T-convergence of the curves of maximal slope for the regularized energies. 

Serfaty’s result on the T-convergence of curves of maximal slope provides a powerful general 
framework. In practice, it can be challenging to verify conditions (i) — (iii). However, for the 
regularized interaction energies E e , these conditions follow from our previous results and the 
following HWI-type inequality, from work by the first author, which allows us to compute the 
strong upper gradients. 


Lemma 5.5 (c.f. 


30 


Suppose pL, v G V 2 ,r( M d ) and E is an interaction energy satisfying 


(E1)-(E2). For ui(x) as in equation (5.2), there exists Cn,d,p,q > 0 so that 

|E (n) - E(u)| ^ 2||VK * g\\ L 2 ifl) d w (g,iy) + C RAP)q ui(d‘( v (g,^)). 


With this, we show that the curves of maximal slope for the regularized energies T-converge 
to the curve of maximal slope of the unregularized energy. 


Theorem 5.6. Suppose the energy E satisfies (E1)-(E2) and the mollifier g> satisfies (Ml)- 
(M3). Let /i e (£): (0,T) -A V 2 , n(M d ) be a curve of maximal slope of E e for the strong upper 
gradient gfig) := 2||Vif £ * h\\l 2 (u) on the metric space (P 2 t R(M d ), dw)- Suppose that /i e (0) is 
well-prepared, in the sense that for some g(0) G V 2 

/i e (0) -A /r(0) weak-* in R(M rf ) and lim E e (/i e (0)) = E(/i(0)). 
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Then for all t E (0,T) ; p £ (t) has a weak-* convergent subsequence p £ {t) ——A p(t), and 
p(t ): (0, T) -a- V 2 ,R.(J& d ) is a curve of maximal slope of E for the strong upper gradient g(p) := 
2\\VK * p\\ L 2 r^. Furthermore, as e ^ 0, 

E e (/x e (t)) -A E(/i(t)) /or all t E [0,T), 

2||VK e * Melina( Me ) -A 2|| VK * mIIl 2 ( m ) in L? oc (0,T), 
and \n' e \ dw -► \p!\d w in L\ oc (0,T). 

Remark 5.7. For any p{ 0) E V2,R(R d ), there exists p, E (0) E V2,r satisfying the conditions of the 
theorem: simply define p e (0) by convolving p(0) with the heat kernel, as in Theorem |3.4[ii). 


Proof. Let ( S e ,d e ) = (S,d) = (V 2 ,R$&- d ), dw)- Note that 

V 2 , R (M d ) = {T£ V 2 (R d ): ||/i|| L oo (Md) sj R} 


is closed with respect to dw, thus V 2 ! n(^k d ) is a complete metric space. Furthermore, any dw 

bounded set of V 2 ,R(J& d ) is relatively compact with respect to weak-* convergence in V(M. d ), 

s 

and its limit points lie in 'P 2 .n{^ d )- Given p e E S 6 , p E S we say p e —^ p if the p e converges 
with respect to the weak-* convergence in V(M. d ). 

We now define the strong upper gradients. For all e > 0, Remark (2.8) ensures that g e (v) = 
2|| VAT e * v\\l 2 (v) is a strong upper gradient of E e on 7^2(M- rf ), hence on V 2 .rt(^k d ). Next, we show 
that g(v) = 2||VAT * ^\\l 2 (u) is a strong upper gradient of E. Throughout, we use the fact that 
if p E L°° (R d ) (1L 1 (R d ), then K *p(x) is continuously differentiable and V(AT * p) = (VAT) *p. 

Suppose u(t) : ( a,b ) -A A^rO^) is an 2-absolutely continuous curve. The function t i-a 
g(i/(t)) is measurable, since it is given by the composition of measurable functions. By Lemma 
[531 


(5.6) | E(u(t)) - E(i/(s))| 2|| VAT * v(t)\\ L 2 ^ t)) dw(v{t), v{s)) + C RAPtq u(dw(v(t), v(s ))). 

we estimate ||ViL * v(t)\\L 2 (v(t)) by breaking VAT into its singular and 


2.8 


As in Remark 


nonsingular parts, 


I VAT * v(t)\\ L 2(v(t)) < ||ViL n * v{t) \\L 2 (y(t)) + ||VAT S * v(t) ||x, 2 (^(i)) 
, 2 \ V2 




VK n (x-y) dv(y,t) 


dv[x,t)\ +VR\\VK 8 * u(t) || L 2 (Hd) 




\ 1/2 

VK n (x - y )| 2 dv{x, t)J dv(y, t) + VR\\VK s \\ L i {Rd) ||p(t)|| L 2 (]Rd) . 


Since X7K n has at most linear growth, the first term is bounded by the second moments of 
v(t), which are uniformly bounded for t E (a, b) since {v(t)}t£( a ,b) is bounded with respect to 
dw- The second term is uniformly bounded since ||i'(£)|| i i = 1 and ||t'(i)||z,°° ^ R- 
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By the absolute continuity of v(t), dw{v{t), v(s)) ^ [' |z/|(r) dr for |z/|(r) E L 2 (a, b ). Com¬ 
bining this with (5.6) ensures that for almost every t 

|E(u(f)) - E(i/(s))| 


E(z/(t)) 


^ lim 

s—yt 


\t — s\ 


/ U m o||V7^ , ..^MI d w (u(t),v(s)) u(d 2 w (u(t),u(s))) 

^ h_™ 2|| VA *v(t)\\ L 2 (u(t)) -- + C R,d,P,q -- 


2||VA^u(t)|| L 2 Mt)) |i/ / |(t), 


where in the last line we use that, for ui(x) defined by (5.2), u(x 2 ) ^ 4xu(x) and lim^-^o uj(x) = 
0. Integrating the above inequality, we conclude, 


|E(u(f)) - E(u(s))| = 


C d 

I ^E(z/(r))dr 


€ 


2||VA * u(r)||x, 2 ( j y(r)) |u / |(r) dr. 


Thus g(v) = 2|| VAT * v^l 2 (v) is a strong upper gradient of E. 

Now we show that for all t £ (0,T), n e (t) -——A fi(t) E ”P 2 ,/i(M rf ), with respect to weak-* 
convergence of probability measures. Since fi € (t ): [0, oo) —> V 2 (J& d ) is a curve of maximal slope 
for the strong upper gradient g e (g) = 2||V-R e * iA\l 2 (ij)-> applying Cauchy’s inequality to (5.4) 
and comparing with (5.5), we obtain (see [l, Remark 1.3.3]) 

(5.7) E e (/j e (s)) - E e (ju e (t)) = J \n' e \ 2 {r)dr 

for all 0 ^ s ^ t < oo. By definition of the metric local slope and the fact that E e (/j, e (t)) ^ 0, 
dw(M0),de(t)) < ( [ \n' e \(r)dr) \g e \ 2 (t) dt ^ T E e (// e (0)) 


for all t G (0,T). Since lim e _;.o E e (// e (0)) = E(/z(0)), the right hand side is uniformly bounded 
for e sufficiently small. Therefore, {/j, £ (t)} e> o is uniformly bounded in V‘ 2 .n(^- d )- and, up to a 


subsequence, fj, e (t) °> fj,(t) G V 2 .n(^ d ) with respect to weak-*convergence in 

It remains to verify criteria (i) — (in) of Theorem (5.3) to conclude that fi(t) is a curve 
of maximal slope of E for g(g,), and the corresponding energies, strong upper gradients and 
metric local slopes converge as e —> 0. 

Criterion (i) follows immediately from Theorem 3.4 part (i) and Remark 3.5 To prove (ii), 
we assume without loss of generality that there exists 0 C < +oo so that 

C = lim inf f \fj/\ 2 (t)dt. 

e -*-° Jo 

Choose a subsequence \jJ' £ \(t) so that lim (E _ s .o fy \£i' e \ 2 (t)dt = C. Then | fi £ \(t) is bounded in 
L 2 (0,s) so, up to a further subsequence, it is weakly convergent to some v(t) G L 2 (0,s). 
Consequently, for any 0 ^ so ^ si ^ s, 


lim [ \p, £ \(t) 

J SQ 


rs i 

dt= v(t) dt. 

j SO 
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By taking limits in the definition of the metric derivative and using the lower semicontinuity 
of dw with respect to weak-* convergence, 

rsi rs 1 

dw(Ve(s 0 ), n e (si)) ^ / 1/41 (t)dt yields d w (//(s 0 ), At(si)) < / v(t)dt. 

J So j So 

By 0 Theorem 1.1.2], this implies that |//|(t) ^ «(£) for a.e. t G (0, s). Thus, by the lower 
semicontinuity of the T 2 (0, s) norm with respect to weak convergence, 
rs ps ps ps 

lim inf / 


e—>0 


/»s /* s p s 

\p' e \ 2 (t) dt = lim / |/4| 2 (t) dt ^ / v(t) 2 dt ^ / |//| 2 (t) dt. 

e ^° Jo Jo Jo 


Finally, we turn to (ra). We assume without loss of generality that 

C = lim inf g e (/j, e ), 

e ->0 

for some 0 ^ C < Too. Choose a subsequence g e {pe) so that lim e _>.o ge(Pe) = C. Since 
lll t e( i )llL°°(M d ) ^ -R for all e > 0 and all t G (0,T), there exists a further subsequence so that 
converges to some limit v in the weak-* topology of L°°(K d ). Since g e also converges to ji 
with respect to the weak-* topology of V(M. d ), for all / G C° 


f(x)v(x)dx = lim / f(x)n € (x)dx= / f(x)p(x) dx\ 


e ->0 


hence, u = /t. 

Theorem 5.4.4] it suffices to show 

lim f f(x)(VK e * He){x) dg e (x) = [ f(x)(VK*g)(x)dg(x) for all / G 
£ ->° Jr^ J]R d 

Then using the convexity and the lower semicontinuity of the function | • | 2 along with the 
weak-* convergence of g e yields the result. 

Following a similar argument as in Lemma |2.5| and defining g e = ip e * * ip € , we have 


; c °°(M d ) 


f(x)(yK e * n e ){x)dn e (x) - / f(x)(VK* g){x) dp{x) 




f(x){VK * p, € )(x) dg e (x) - / f{x){VK * g)(x) dg(x) 


f(x)[(X7K * fl e )(x) - (V/\ * y)(x)] dg e (x) 


+ 


f{x){VK * n){x) dg t (x) - / f(x)(\7K*g)(x)dg(x) 


— : A € + B e 


Since —> p weak-* in V(M. d ) and f(x)(VK * p)(x) is bounded and continuous, lim e _>o B e = 0 

It remains to show lim e _>o A e = 0. First, note that p e — > p in the weak-* topology of L 
as e —> 0. Indeed, for any / G L 1 (M d ), we have that 

/ f(x)dp e (x)~ / f(x)dp(x) = / ip e * f * <p e (x) dp e (x) - / f(x)dp(x) 

■J R d J R d JR d J R d 


^ R\\ip e * f*<p e - f\\ L i + 


f(x) dp e (x) - / / (x) dp(x) 
















REGULARIZED NONLOCAL INTERACTION ENERGIES 


25 


where both terms approach zero as e —> 0. Returning to A e , 

A < [ \f(x)(VK * (jl e -fl))(x)\ \/j e (x)\ dx 

JR d 

<-R||/IIl°°(r<*) / A R * (A - n)(x)\dx. 

J supp / 

Since the integrand has at most linear growth, it is bounded on the compact set supp /, and we 
may apply the dominated convergence theorem, provided the integrand converges pointwise. 

When 0<g^l,VRTis the sum of a bounded continuous function and an integrable 
function, and since °> fl in both weak-* probability and weak-* L°°(R rf ), the integrand 

converges for each x. On the other hand, when 1 < q ^ 2, it suffices to show that X7K is 
the sum of a continuous function, which is uniformly integrable with respect to jl e , and an 
integrable function [lj Lemma 5.1.7]. In particular, is enough to show that |x| 9_1 is uniformly 
integrable with respect to fl, e . Since \x\ ^ k implies that \x\/k ^ 1 we have that 

lim f |a?] 9 " 1 d/2 e ^ lim — f \x\ g djl e ^ lim — E(/2 e ) = lim —E e (// e ). 

k—yoo Jjxj^k k—too k J\ x \^k fc—>• oo k k —»oo k 

Since fl £ (t) is a curve of maximal slope for E e , ( |5.7| ) gives E e (/j, e (t)) ^ E e (// e (0)) for all t G (0, T). 
The well-preparedness of the initial data gives E e (/r e (0)) < C, hence the result follows. □ 


Remark 5.8. The extension of Theorem |5.6| to energies defined via kernels satisfying the more 
general conditions (H1)-(H5) from section [4] requires also a generalization of the technical 
result of Lemma 5.5, which plays a crucial role in the proof of the convergence theorem. As 


such a generalization is beyond the scope of our current study we leave this to future studies. 
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